SPhT-T06/012 



Thermal flow in the gravitational 0(n) model 



Ivan K. Kostov 12 

Service de Physique Theorique, CNRS - URA 2306 
C.E.A. - Saclay, 
F-91191 Gif-sur-Yvette, France 



Abstract 

We study the massless flow from the critical point (dilute loops) to the low-temperature 
phase (dense loops) of the 0(n) loop gas model when the model is coupled to 2D gravity. 
The flow is generated by the gravitationally dressed thermal operator $1,3 coupled to the 
renormalized loop tension A ~ T — T c . We find that the susceptibility as a function of 
the thermal coupling A and the cosmological constant fi satisfies a simple transcendental 
equation. 
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1 Introduction and summary 



It is well known that the critical phenomena on flat and fluctuating lattices are deeply related. 
For each critical point described by a "matter" CFT, the "coupling to gravity" consists in adding 
a Liouville and ghost sectors and dressing the scaling operators by exponents of the Liouville 
field. In such simplest theories of 2D gravity are quite well understood nowadays due to the 
progress towards the exact solution of Liouville theory achieved in the last decade. 

On the other hand, almost nothing is known about theories of 2D gravity in which the matter 
field has massless excitations but is not conformal invariant. A typical example of such a theory 
is the 0(n) loop model on a honeycomb lattice [1]. This model has two non-trivial critical points, 
the dense and the dilute phases of the loop gas, described by two different CFT's. The massless 
flow relating these two CFT's is generated by the thermal operator $1,3. Other theories with 
massless flows have been studied in [2-4]. 

The thermal flow in the gravitational 0(n) model is expected to be described by Liouville 
and matter CFT's coupled through the operator gravitationally dressed thermal operator $1,3. 
It is not known at present how to solve such a theory. On the other hand, the corresponding 
microscopic theory, the 0(n) model on random planar graphs, can be solved exactly using its 
dual formulation as a matrix model [5-8]. The exact solution in the case of general potential 
has been formulated in [9,10]. 

The aim of this paper is to work out, using the matrix model formulation, explicit expressions 
for the partition functions of the gravitational 0(n) model on the disk and on the sphere along 
the thermal flow. In this section we give a brief introduction in the 0(n) model and present our 
results. 

The 0(n) model can be defined on any trivalent graph Q. The local fluctuating variable 
associated with the sites r E Q is an 0(n) spin with components Si(r), . . . S n (r), normalised so 
that Tr S a (r)Sb(r') = 5 a b5 rr >. The partition function on the graph Q is defined as 

Z 0(n) (G;T) =Tr [] (l + ± £ S a (r)S a (r')) (1.1) 

<rr'> a 

where T is the temperature and the product runs over all links <rr'> of the graph Q. Expanding 
the trace as a sum of monomials, the partition function can be written as a sum over all 
configurations of self-avoiding, mutually avoiding loops that can be drawn on Q, 

Z o {n) (9, T) = J2 T ~ Ctot nM °° ps > (1-2) 

loops 

as shown in Fig.l. Here £ t ot is the total length of the loops, equal to the number of occupied 
lattice edges and M\ oops is the number of loops. Unlike the original formulation, the loop gas 
representation (|1.2j) makes sense also for non-integer n and has a continuous transition in two 
dimensions for \n\ < 2. In this interval the number of flavors can be parametrized as 
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Fig. 1: Loops on a trivalent fat planar graph 



n = 2cos7ri/, < v < 1. (1.3) 

The phase diagram of the loop gas on the infinite regular trivalent graph, the honeycomb lattice, 
was first established in [1]. At the critical temperature T c = 2cos^ the loop gas model is 
solvable and is described by a CFT with central charge 

c = 1-6^-. (1.4) 

critical 1 + V ^ ' 

For T > T c the theory has a mass gap. The low-temperature, or "dense", phase T < T c is a flow 
to an attractive fixed point [18] at T^ ense = 2sin?f, where the theory is again solvable and is 
described by a CFT with smaller central charge 



c. =1-6^. (1.5) 

dense 1 — t/ V / 

The scaling behaviour of the model in the vicinity of a critical point is described by an action 
of the form [4] 

A = Critical + 5T [ $ 1)3 (1.6) 



where 5T = T — T c and $13 is the thermal operator with conformal dimensions 

Ai )3 = A 1>3 = (1.7) 

Microscopically the thermal operator $1^ counts the total length of the loops £tot- Added to 
the action, it generates a mass of the loops. For ST > the deformation Ifl.fiJ) describes, going 
from short to the long distance scales, the flow to a massive theory with mass gap 



m ~ 



ST l/p-2Ai, 3 ) =ST *&_ M 



When 5T < the deformation (jl,fi[) describes a massless flow between two different CFT with 
central charges (|1.4|) and ()1.5|) . In the CFT for the dense phase, describing the IR limit, the 
flow to the attractive critical point is generated by another, irrelevant, operator, ^3,1- 

The 0(n) model on a fluctuating lattice [5] can be formulated as a statistical ensemble of 
trivalent planar graphs covered by self- avoiding and mutually avoiding loops. The critical ther- 
modynamics is controlled by the temperatore T and one extra parameter, the bare cosmological 
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constant k coupled to the size of the graph. For example, the partition function on the sphere 
is defined as 



^K,T) = ^K A Z 0{n) (g,T) (1.9) 
Q 

where the summation is taken over all connected fat graphs Q with the topology of a sphere and 
A is the area (the number of vertices) of the graph. 

As in the case of a flat lattice, the model has three critical points. Each critical point is 
characterized with the "string susceptibility" exponent 7 str related to the matter central charge 
by c mat ter = 1 — 67g tr /(l — 7 s tr) • [H]- Qualitatively, at the critical point k = k*,T = T* both the 
area of the graph and the length of the loops diverge. The flow to the massive, high-temperature, 
phase is along the critical line k = ki(T) where the area of the graph diverges, while the loops 
remain finite. The flow to the dense, low-temperature, phase is along another line, k = ku(T), 
where the area of the graph diverges because of the diverging length of the densely packed 
loops. Thus the continuum limit is described by the vicinity of the critical line k = k c (T), which 
consists of two branches meeting at the critical point: 



K = K C (T) 



kj(T), T > T*; 
kh(T), T < T*. 



The exponent 7 str is given in the three phases by 



7str 



( 1 

2 



__ v 
\ l-u 



high-temperature phase (massive loops); 

critical point (dilute loops); (1-10) 
low-temperature phase (dense loops). 



In the continuum limit the loop gas model on a fluctuating lattice is described in terms of the 
renormalized coupling constants 

fi~K c {T)-K 2 , A ~T*-T. (1.11) 

Denote by <&(/U, A, hb) the partition function on the disk with boundary cosmological constant 
Hb- It is given by the Laplace transform of the disk partition function <E>(/x, X,£) with fixed 
boundary length I: 

POO 

$(Ai,A,Ai B )= / d£e-^ e $([i,\,e). (1.12) 
Jo 

An important characteristics of the model is the boundary entropy — M, defined as the 
exponent in the exponential decay of the disk partition function when the length tends to 
infinity: 

M = — lim kg (1.13) 

As each loop can be considered as two boundaries glued together, the effective loop tension is 
equal to twise the boundary entropy. 

Summary of the results: 
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We have found that the susceptibility x ~ is related to the boundary entropy as 

X = M 2u . (1.14) 



The function M = M(A, (i) is determined by comparing the expressions for the derivative d^l^g 
and d^g&n, which we obtained by solving the saddle point equations for the 0(n) matrix model. 
The derivatives of the disk partition function are given in parametric form by 

[is = Mcoshr, 

fy B $U = ~ 2Ml+V cosh (l + v)t 

+2AAf 1 " iy cosh(l - v)t, 

9 m *| Mb = v~ x M v cosh vt. (1.15) 

(The normalization of the coupling constants is chosen for our convenience.) The compatibility 
of these two expressions implies the following transcendental equation for the boundary entropy: 

\k = (1 + u)M 2 + AM 2-2 ". (1.16) 

Equation Q1.16JI generalizes the previously obtained results for the scaling behavior of the bound- 

i i 
ary entropy in the dilute and the dense phases, M ~ fx? in the dilute phase and M ~ A 2 - 2 " in 

the dense phase [12]. 

As the susceptibility is related to the boundary entropy by (|1.14|) . it satisfies the transcen- 
dental equation 

li={l + v) x « +A X 1 ^ J (1.17) 

which coincides, after a redefinition of the variables, with the eqation found in [13] for the 
susceptibility of the gravitational sine-Gordon model. 

The dimension of the coupling A ~ fj> u matches the gravitational dimension 613 = 1 — 
v, obtained from Aj.,3 by the KPZ scaling relation A = S ^t^ [11]- This is consistent with 
the conjecture that for finite \/[i u the theory is described by a perturbation of the critical 
point A = by the Liouville-dressed thermal operator $1,3. An analysis of eq. (|1.17|) in the 
framework of Liouville gravity was performed by Al. Zamolodchikov [14]. After integration 
one can reconstruct the expansion of the partition function T in A, which generates the n- 
point correlation functions of the thermal operator. The result for n < 4 matches with the 
calculations recently published by A. Belavin and Al. Zamolodchikov [15] as well as with the 
formula conjectured in [16] on the basis of the ground ring identities. We leave the analysis 
of the case of a strong perturbation from the perspective of the worldsheet CFT to a future 
publication [17]. 

2 Solution of the 0{n) matrix model in the planar limit 

2.1 The 0(n) matrix model 

The partition function (|1,9|) of the 0(n) model on random triangulations is generated as the 
perturbative (t'Hooft) expansion of the free energy of a model of n + 1 hermitian N x N matrix 
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variables X and Y = {Y\, Y n } [5]: 

>Jdxf[dY a e NK ~ 2 (-^ ^-|E: = x^+^ 3 +^E: =1 ^ 2 ). (2.18) 



a=l 



We are interested in the planar limit N — > oo, where only genus zero planar Feynman graphs 
survive. Each such Feynman grah realizes a loop configuration in the sum l|1.9[) . 

Integrating out the K- variables and shifting X — > X + one obtain one-matrix integral of 
the form 

= J dX e~P tr v{x) |Det (1 ® X + X l)\~ n/2 (2.19) 

where f3 = N/k 2 and the coefficients of the cubic potential 

V(X) = hX + t 2 X 2 + t 3 X 3 (2.20) 
are expressed in terms of the temperature T as 

h = ^n, * 2 = ^, t 3 = -i (2.2i) 

It will be convenient to take the large N limit by sending (3 to infinity keeping the ratio k 2 = f3/N 
finite. 

After diagonalization gXg~ l = {x\, ...,xjs?}, g G SU(N), the matrix integral (|2.19j) can be 
reformulated as a two-dimensional Coulomb gas of N charges constrained at the real axis: 



e 



, N N N 

T ~ / TT dxi e-^^ JJixi - x 3 ) 2 TT ( Xi + Xj )S. (2.22) 



i=l i<j M=l 



In the limit N — > oo the integral is saturated by the saddle point described by the classical 
spectral density p(x) supported by the interval [—a, —6] on the negative semi-axis, where < 
b < a. We will normalize the spectral density as 

r b 

I dxp(x) = ^ = K 2 . (2.23) 

J — a 

The spectral density is determined from the saddle point equation 

-ft -ft 

V'(x) = 2V [ dy-^- - n [ dy-^-, x G [-a, -61. (2.24) 
J x -y J x+y 

—a —a 

The edges of the eigenvalue distribution are functions of the external potential, i.e. of the 
temperature T and the cosmological constant k. 

The disk partition function with boundary length I is defined as 

/•-ft 

= Jtr e £X = \ / e M p(\)d£. (2.25) 

J — a 

The position of the right endpoint, —6, determines the large t asymptotics of the disk partition 
function, $>(£) ~ e~ M , and has the statistical meaning of boundary entropy for unit length. The 
boundary entropy is negative because shifting the matrix variable X we effectively performed a 
subtraction. 
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Knowing the position of the edge b = b(T,K), the equations of the two branches of the 
critical line mentioned in the Introduction are obtained as follows [5]. The branch n = k,jj(T) 
is determined by the condition that the right endpoint reaches the origin: 

6 = => k c = kjj(T). (2.26) 

When b > the loops are not critical, the critical behavior is that of pure gravity (c mat ter = 0) 
and the position d K b ~ d\T must behave as (k c — k) -1 / 2 . The branch k c = k/(T) is determined 
by the condition that the position of the endpoint b develops a singularity in k: 

db 

— -^oo => k c = Ki (T). (2.27) 
Ok t 

Finally, the critical point {T, n} = {T*,k*} is the common endpoint of the two critical lines: 
k* = kj(T*) = k h (T*). 



2.2 Functional equation for the loop field 

The saddle point equation can be solved directly, see e.g. [6] for the solution in the special case 
6 = 0. It is however more advantageous to reformulate the problem in terms of the collective 
field, or loop operator, 

N 

Hx) = -^log(x- Xi ), (2.28) 
i=l 

which is the Laplace transform of the disk partition function (|2.25[) . The geometrical meaning 
of the operator $ is that it creates a boundary with, in general complex, boundary cosmological 
constant \i% = x. The vacuum expectation value ( $(x) ) is equal to the disk partition function, 
the connected correlator ( <&(x)&(x') ) gives the partition function on a cylinder, etc. 

The saddle point equation Q2.24JI can be reformulated as a boundary condition for the current 

J(») = - 2V ' ( 'g H -W * GC ' ( 2 - 29 ) 
which is analytic on the complex plane cut along the interval [—a,— b}. Namely, the values of 
J(x) on both sides of the cut are related by 

J(x -HO) + J(x - i0) +nJ(-x) = xe[-a,-b]. (2.30) 

The function J(x) is completely determined by the boundary condition (|2.3Uf) and the first four 
coefficients of its Laurent expansion at infinity, 

3 

J(x)= ^ Jnoc'" 1 " 1 = J 3 x 2 + J 2 x + Ji + Jox- 1 + [J(x)] <0 , (2.31) 

71= — oo 

which are expressed in terms of T and k, 

J = k 2 , Ji = - T $,Zn)> J 2 = ^3 = 2^- (2.32) 
We will split, as in [8], the function J(x) in two 'chiral' pieces, 

J{x) = J + {x) + J_(x), (2.33) 
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satisfying simpler boundary conditions along the cut: 

J+(-*±*0) = -e±™J + (x±iO) (b<x<a) (2 34) 

J_(-s±*0) =-e^ mv J-{x±iO) ^< x < a >- 

The phase v is defined by The boundary condition 1)2 .3U|) . respectively 1)2.34)) . implies a 

quadratic functional identity for the current [12] 

[J\x) + J 2 {-x) + nJ{x)J{-x)] <Q = (2-n 2 )[J+(x)J_(x)] <0 = 0, (2.35) 

where [ ]<o denotes the negative piece of Laurent expansion. Therefore J+(x)J-(x) is an even 
polynomial of 4 degree whose coefficients are known functions of k and T. In order to avoid 
heavy expressions we rescale x and b as 

x = x /a, b = b/Aa (2.36) 

and use directly the Laurent expansion of the chiral components J±: 

J ± ( X ) = e TW2 (j- h k x 2k ~ l ± i 'hk+i xA . (2.37) 

V k k J 

The new expansion coefficients are related to the old ones by 

^ 2k = " s /2-n ^ 2k ' ^2k+l = J 2+n J2k+l- (2.38) 

Then (|2.35|) implies the functional equation 

J+(x)J_(x) = A + Bx 2 + Cx 4 (2.39) 

where the A, B and C are expressed in terms of J3, J2, Jo and the first moment of the eigenvalue 
density W\ = J_i: 

A = 2W 1 J 3 + Jl + 2J J 2 , 5 = 2JiJ 3 + j|, C = Jl (2.40) 
2.3 A criterium for criticality 

The functional equation 1)2. 39|) leads to an algebraic equation for J(0): 



J 2 (0) = (2-n) J+(0)J_(0) = (2 - n) A => J(0) = ^(2 - n)A (2.41) 

The value J(0) is analytic function of the couplings k 2 and T in the domain where the loop 
gas partition function is convergent. The critical lines are therefore given by the border of the 
domain of analyticity of J(0). The latter is given by the condition A = 0, or 

21¥iJ 3 + Jl + 2J0J2 = 0. (2.42) 
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2.4 Solution along the critical line k = Kn(T) 

The Riemann surface of J(x) consists of infinitely many sheets, except for the case when v is a 
rational number. The first, physical, sheet has one cut [—a, —b] while all the other sheets have 
two cuts [—a, —b] and [b, a] (Fig. 2). If we find a global parametrization of the Riemann surface 
that resolves all branch points, then the boundary conditions (|2.^4j) will become quasi-periodicity 
conditions and can be solved (in our case) in terms of theta functions. 



P(x) 



Fig. 2: The Riemann surface of J(x) 

To illustrate the method, let us first reproduce the solution along the critical line T = T c (k), 
found originally in [6] by applying the Wiener- Hopf method to the integral equation (|2.24j) . 
Along this critical line b = and the infinity of simple ramification points at x = ±6 merge 
into a single ramification point of infinite order at x = 0. The Riemann surface is then globally 
parametrized by a hyperbolic map 

^ - dh7 (2 - 43) 



where the points with parameters s and —s must be identified. The boundary conditions (|2.34jl 
become a quasi-periodicity conditions in the complex s-plane: 



J+(s ± iir) 
J_(s ± iir) 



-e^™ J_(s). 

The unique solution of (|2.44|) with the asymptotics (|2,37|) is 



J±(s) = di 



cosh s 



+ d 2 



e ±(u-l)s 

cosh 2 s 



with 



d2 = J3, di = J 2 - (1 - v) J 3 . 
One can check that the quadratic equation (|2.39l) is satisfied with 



(2.44) 

(2.45) 
(2.46) 



,4 = 0, B = d 1 {2d 2 + di), C = d\. 
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Comparing with the asymptotics (j2.37j) we find two more relations 



Ji = yh-^h, Jo = \J2- 1 =^-h, (2.47) 

which determine the left branchpoint a = a(«) and the critical line k c = kh(T) that describes 
the dense phase. The line k c = kh(T) starts at the critical point T* where d\ = 0. 

In a similar fashion we calculate the Jo-derivative of J(z), which we rescale for convenience 

as 

d J{x) = aV2-n d J(x). (2.48) 

The meromorphic functions d$J±{x) satisfy the same boundary condition Q2.56JI but have less 
singular behavior at infinity: 

d J ± {x) = ±ie^l 2 ( K l^i^ + ...), x^oo. (2.49) 

The only solution of 1)2.56(1 with such asymptotics is 

d J± = ±ie TUS coth s. (2.50) 



We finally obtain for the current (|2.29|) and its derivative 

J = 2d!^ + 2(i 2 cosh(1 -^ 



coshs 2 cosh^ , 

BqJ = sinhz/s coths. (2-51) 



The scaling behavior for \x\ <C 1 is 



J{x) ~ (p-b) 1+u - A (p,b) 1 ~ p , 

d J(x) ~ {jinf. (2.52) 



where 

\ — hzR -h. 



2.5 The general solution 



A solution of the 0(n) model in the most general case is presented in [9, 10]. However, we found 
easier to perform an independent computation rather than to use the results of [9, 10]. 

For generic values of the couplings there exists a global parametrisation in terms of Jacoby 
elliptic functions. The uniformization map u — * x(u) is defined as 



dn(n, k) ' 



k' = VI - k 2 = i (2.54) 



(our notations are those of Gradshtein and Ryzhik [20]). 

The map has the symmetries x(u) = x(—u) = x(u + 2K) = —x(u + 2iK') and parametrizes 
the physical sheet by the rectangle [0,K] x [—2iK',2iK'] in the w-plane. The whole Riemann 
surface is parametrized by the orbifold of the -u-plane with respect to the symmetry u — > —u. 

The functions J±(u) must be periodic in u — > u + 2K 

J±{u + 2K) = J±{u) (2.55) 
and, by the boundary condition (|2.34j) . quasi-periodic in u — > u ± 2iK': 

J+(u ± 2iK') = - e ±in,/ J + (u), 
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J_(u±2iK') =-e^ v J + {u). (2.56) 

The general solution of these equations is given in terms of the theta function Y w {u) defined as 
follows: 

The map (|2.54j) is a particular case of this function, x(u) = bY\{u). We will need the following 
properties of the function Y w (u): 

1) Symmetries and (quasi) periodicity: 

YUu) = y_ w (-«) = Y„ +2 {u) = YUu + 2K) 

Y u {u ± 2iK') = e ±i7Tl ^Y u (u) (2.58) 

2) Quadratic relations 

y w (u)y_ w (w) = 1 + sd 2 (i/K) (x 2 - k' 2 ) , (2.59) 

yi_ w («)yo,(«) + y_i4*,(«)y-«(«) = £ £ (2.60) 

yi- w («)y w («) - k-i-^M^M = f^MO ft,*, (2.61) 

3) Particular values: 



u 




y w 





6 


1 




a 


nd(wK) 


K±iK' 


00 


00 


K ± 21KT' 


—a 


e ±im "nd(u;if) 


±2iK' 


-6 





4) Asymptotics at x — > oo: 

r„(z) = sdMQ (x - + ...) (2.62) 

One can easily check that the solution of (|2.56|) in the case of cubic potential is given by the 
linear combination 3 

J±{u) =D l Y ±(l ^ ) {u) + D 2 Y TV {u)x{u). (2.63) 

Since J+(— u) = J-(u), the function J = J + + J_ is even function of u, hence a function on the 
orbifold. The coefficients D\ and D 2 and the branch points a and b = k'a are obtained by from 
the asymptotics at infinity given by (|2.31[) . To simplify the expressions, we will introduce the 
rescaled coefficients 

d l = D 1 ^> , d 2 = D 2 sd(uK) . (2.64) 

3 The solution for a generic polynomial potential is obtained by replacing the coefficients Di and D2 with entire 
functions of x 2 [9], 
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Comparing the coefficients of the two leading of the expansion at x — > oo we find 

d x = J 2 - f(vK)j 3 , d 2 = J 3 . (2.65) 

Note that (|2.46j) = (|2.65j) \k'=o- We need two more relations that determine the positions a and 
b of the branch points. Instead of expanding further around the point x = oo, it is simpler to 
use the fact that the solution (|2.6Hj) satisfy a quadratic relation whose l.h.s. is identical to that 
of (|2.39[) . This quadratic relation follows from (|2.61|) and (|2.6Uj) : 

J+ J_ = d 2 x 2 (x 2 + cs 2 {uK)) + 2d 1 d 2 ^(vK) x 2 

+d 2 (x 2 +k' 2 sc 2 {vK)). (2.66) 

Taking x = we obtain for the coefficient A 

A = 2WiJ 3 + Jl + 2J J 2 = k' 2 sc 2 {uK)d\. (2.67) 

The condition for criticality A = is achieved either if k' = (dense phase), or if d\ = (pure 
gravity), or if k' = d± = (critical point, or dilute phase). To evaluate the boundary entropy 
M ~ b(X, /«) we will proceed as follows. We will first calculate independently J and d^J, then 
take the continuum limit and compare the second quantity with the derivative of the first in fi. 
This will give an equation for the function b(\,fi). 

In a similar way we evaluate the derivative dJ±(x) as the unique solution of (|2.56f) and 



2.49j) . It is given by 



where 



§ J± = ±ids(uK) (2.68) 
d u x 



d u x = i\J (x 2 — l)(x 2 — k' 2 ). 
The solution (|2,68|) leads to a simple formula for dod u &: 

dod u $ = —6s(uK) (Y v — Y- v ). (2.69) 

Finally, comparing the subleading terms of the expansions (|2.49|) and (|2.62j) , we find the deriva- 
tive of W\ = ^(tr M), 

whose singular part gives, up to a normalization, the susceptibility \ ~ —d 2 T. 
2.6 Continuum limit 

The continuum limit a — > oo is achieved when both x = x/a and b = b/4a = A//4 are small. In 
this limit the uniformization map (|2.54|) degenerates to 

x(u) = bcoshu (2-71) 

and the theta function ()2.57|) is approximated by 

Y±Uu) = e^ u . (2.72) 
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The two coefficients of the solution (|2.63|) are now given by 

D 1 = 2b 1 - u d 1 , D 2 = 2b u d 2 (2.73) 

and for the current J = J+ + J_ we find 

J =4J 3 (-A6 1 - !y cosh(l-z^)u + 6 1 +^cosh(l + i/)n) (2.74) 

with the coupling A defined earlier in 1)2.53)1 . On the other hand, from 1)2.68)1 we find for the 
derivative do J 

d J = V- 1 (2.75) 
Now let us compare 1)2.75)1 with the derivative of 1)2.74)) : 

d J\ x = d b (J - j^ll) • J 

= 4J 3 Bob (-\(l-v)b-"-(l + v)V') (2.76) 

V / smn u 

The two expressions coincide under the condition 

Jo - J c o = -2 J 3 (\b 2 -^ + (1 + v)b 2 ) , (2.77) 



where the T-dependent integration constant Jq corresponds to the value k c = kjj(T). 

Eq. ()2.77)l is identical to the transcendental equation ()1.16)l with the cosmological constant 
H and the boundary entropy — M being normalized as 

H = M = b. (2.78) 

it/3 

The susceptibility \ = —d 2 J- is given, up to a constant factor, by the leading singular term in 
the small k' expansion of ()2.7U)) : 

-doW 1 = l-v + 2b 2u + ... (2.79) 

Since we have not yet normalized the string interaction constant g s ~ j3, we have the freedom 
to normalize the free energy. We do it so that \ = M 2l/ . 
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